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Abstract 

We describe non-conventional localization of the midband E = state in 
square and cubic finite bipartite lattices with off-diagonal disorder by solving 
numerically the linear equations for the corresponding amplitudes. This state 
is shown to display multifractal fluctuations, having many sparse peaks, and 
by scaling the participation ratio we obtain its disorder-dependent fractal 
dimension D2. A logarithmic average correlation function grows as g(r) ~ 
7/lnr at distance r from the maximum amplitude and is consistent with a 
typical overall power-law decay |V>( r )l ~ r ~ ,? where rj is proportional to the 
strength of off-diagonal disorder. 
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It is well-known that off-diagonal disorder in the nearest neighbors of one-dimensional 
(ID) lattice is responsible for anomalous localization of the corresponding state at the band 
center E = [|l]-|4j]. The log-amplitude In a ^ site % executes a random- walk and the typical 
wave function decays exponentially at distance r from its maximum as |^(r)| ~ exp(— a-y/r), 
with a a disorder-dependent parameter. The E = is a special state arising from a sub- 
lattice symmetry, which can exist in the presence of off-diagonal disorder and causes the 
Hamiltonian to change sign under the transformation — > (—1)*^. This symmetry exists 
in bipartite lattices with two sublattices, one connected to the other 0,0], and is also known 
as chiral symmetry. The corresponding eigenvalues appear in pairs with energies (E, —E) 
and since the band center E = is a special energy the corresponding wave function can be 
easily constructed. 

In 2D and 3D bipartite systems, such as square or cubic lattices, the question of local- 
ization at the band center with off-diagonal disorder is not resolved as in ID and the decay 
of the wave function amplitude from its maximum is not yet firmly established. It is be- 
lieved that non-exponential decay occurs in 2D and even absence of decay in 3D where 
extended states exist for weak disorder. The approach of predicted exp(— jyhxr) decay 
in every dimension, which is rather weaker than a power-law. The theoretical interest in 
this problem has been recently revived. On one hand, due to the related problem of random 
Dirac fermions [0,[8| which, however, starts from a different zero disorder limit. In the case 
of off-diagonal disorder one has a Fermi surface with many points forming a square while 
in the problem of Dirac fermions without disorder at zero energy four points exist only. In 
the former case the pure density of states at E — has the well-known two-dimensional 
log-singularity while in the Dirac case the density of states approaches zero at the band 
center ||. On the other hand, the study of localization with off-diagonal disorder could 
be important for understanding properties of several realistic 2D systems, such as states 
in semiconductor quantum wells and high mobility silicon MOSFETs pf, current-currying 



states in quantum Hall effect [[TIJ, etc. 

In this Letter we present an exclusive study of the midband E = state in 2D (squared) 
and 3D (cubic) bipartite lattices with disorder in the nearest-neighbor hoppings. We com- 
pute the coefficients of the wave function for large finite systems and obtain the participation 
ratio. The amplitudes are shown to display multifractal fluctuations with many scattered 
peaks, characterized by a fractal dimension which depends on the strength of disorder. An 
appropriate correlation function is shown to behave as g(r) ~ rj In r where rj also depends on 
the strength of disorder. It implies a power-law decay |V ; ( r )| ~ r ~ v of the peak heights from 
the maximum peak, in contrast to the one-dimensional typical decay |'0( r )| ~ exp(— «i/r). 

We consider a tight-binding Hamiltonian in 2D, 3D lattices with random nearest-neighbor 
hoppings 

ff = £(Vta + H.c.), (1) 

(ii) 

where q is the annihilation operator of an electron on the ith site and the nearest-neighbor 
hopping integrals ty- are real random variables satisfying the probability distribution 

1 W W 

PQnUj) = — for - y < lnfy < —. (2) 
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This disorder, commonly used in ID, guarantees positive hoppings and the typical variation 
of the random elements W is believed to be a good measure for the strength of off-diagonal 



disorder in any dimension JTI]]I2]] . If, instead, we choose to distribute directly between 
— W/2 to +W/2 the disorder can never be made strong enough since there is no energy scale 
in the Hamiltonian. In the following we consider off-diagonal disorder form Eq. (2) and 
ignore diagonal disorder by setting the site energies equal to zero. 

A bipartite lattice consists of two interconnected sublattices A and B with the random 
hoppings tij connecting the sites of one sublattice to the sites of the other. It was shown |J 
that a finite system with ua sites on sublattice A and n B sites on sublattice B {n B > ua) 
has at least n B — Ua linearly independent eigenfunctions with eigenvalues exactly E — 0. 
Moreover, the amplitudes of these states vanish on the sites of sublattice A. In the present 
study we consider finite square and cubic lattices with their L-site edges along the main 
directions in 2D and 3D. In the calculations L is odd and Dirichlet boundary conditions are 
applied with sublattice A having (^^-) sites and sublattice B has (-^4^-) sites 
for the square (cube). According to Ref. U only one E = state exists in this case with 
finite amplitude on sublattice B and zero amplitude on sublattice A and the corresponding 
amplitudes satisfy the linear system of equations 

tjj+sipj+s = 0, for all j G A, (3) 

8,(j+SeB) 

where 5 is summed over the nearest neighbors of site j so that j + 6 belongs to the sublattice 
B where the amplitude is finite. The total number of equations Z/2 ~ 1 { LA 2 l ) derived from the 
A sites allows to obtain the {^p-) amplitudes on the B sites. This procedure suffices 
to determine uniquely the E = state since the the normalization condition J2i \ipi\ 2 — 1 
accounts for the one missing equation. 
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Fig. 1. The logarithmic amplitude In |-0 1 of the E — wave function in 2D squared lattices with 
off-diagonal disorder strength (a) W = 0.1, (b) W = 1 and (c) W = 10. The darkest regions cover areas 
with the higher In \ip\, between and —1, and the lightest regions areas with the lowest In \ip\, between —12 
and —13, while the various degrees of gray denote 13 linear scales. 

We have solved numerically the linear system of equations (3) for different strengths W of 
the off-diagonal disorder. In Fig. 1 typical pictures of the logarithmic amplitudes of the E = 
wave function in 2D are shown which display fractal characteristics. Similar fractal patterns 
are obtained for all examined disorder strengths in 2D and 3D. The logarithmic amplitude for 
very weak disorder (small W) is periodically distributed in Fig. 1(a), the fractal character 
of the state is obvious for intermediate disorder W in Fig. 1(b) and for higher values of 
disorder (large W) the area with significant amplitude becomes vanishingly small fraction of 
the total in Fig. 1(c). The latter case implies stronger decay of the E = state, similarly to 
what one expects for ordinary localized states with strong diagonal disorder. However, the 
important difference to conventional localization is that the maximum is not concentrated 
in a small region of space, but many maxima of almost similar heights exist. These peaks 
with relatively large amplitude are randomly scattered over space and exist also for large 
disorder (Fig. 1(c)). Therefore, the E = state in the presence of off-diagonal disorder 
shows unusual localization properties. 

Firstly, we investigate the fractal structure of the normalized wave function by computing 
the inverse participation ratio 

P = Y.m'-L- D \ (4) 
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where the scaling of p with size L defines the fractal dimension D 2 . It is well-known that 
extended states are space filling with D 2 = d, the space dimension, and localized states are 
point-like with D 2 = 0. The probability distributions of hip for various finite-size L systems 
over a sufficient number of disorder configurations are shown in Fig. 2 to be similar to the 
ones reported in ||13|| . For example, in Fig. 2(b) we plot the distributions for the intermediate 
disorder strength W — 1, which are roughly invariant in shape as the size L gets larger and a 
sudden drop occurs for small values of hip with a tail for large In p. We focused on the peaks 
of such distributions by determining p max from the maximum probability of p as suggested 
in [13]. The p ma x shift towards lower values by increasing the size so that the variation of 
this typical value of p according to Eq. (4) is shown in the insets of Fig. 2 where we plot 
ln(p max ) vs. InL. The fractal dimension for W = 1 is D 2 ~ 1.55. 

The behavior of this special state for off-diagonal disorder is, somehow, reminiscent of 
the conclusions reached in [14J] where the 2D wave functions are multifractal in the presence 
of disorder for length scales shorter than the localization length. For off-diagonal disorder 
a similar situation is encountered since the usual (exponential) localization length is guar- 
anteed to diverge at E = 0. In Fig. 2 we plot the probability distribution of the inverse 
participation ratio for W = 0.1, W = 1.0 and W = 10, where D 2 is well-defined from the 
slope of the \np max vs. InL almost linear curve in the investigated sizes. For small disorder 
W = 0.1 we find D 2 pa 2, which implies an almost extended state, on the contrary, for 
strong disorder W = 10 the logarithm of the typical inverse participation ratio is close to 
zero, which implies strong localization with D 2 pa 0. 
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Fig. 2. The probability distribution for the logarithm of the inverse participation ratio for the E = 
wave function in square (2D) lattices with sizes L = 23; 35; 71; 107 and 121. The inset shows \n(p max ) vs. 
InL for (a) W = 0.1, (b) W = 1, (c) W = 10. 

Finally, we examine the typical decay of the E = state by introducing an appropriate 
correlation function. Following the definition of we have studied 

^(r) = (ln^J),r=|r,-r w J, (5) 

where i max is the site of the lattice with the largest amplitude and the average (...) is 
taken over a sufficient number of random configurations and various spatial directions. In 
ID bipartite lattice, where the log-amplitude of the E = state with off-diagonal disorder 
exhibits random- walk behavior, the correlation function is g(r) ~ a/t for large r |. In 
Fig. 3 we plot g(r) vs. lnr in 2D by increasing W. These curves for not too small r can 
be reasonably well fitted to the relation g(r) = rjlnr + (3, with parameters rj and (3. For 
example, for W = 1 we find rj = 0.44 ± 0.01, f3 = 1.51 ± 0.02 and the average amplitude 
asymptotically decays from its maximum as a power law |"0( r )| ~ r ~ v - We conclude that the 
multifractal amplitude pattern shows many random peaks and the power-law decay reflects 
the long-range decay of the peak heights from the highest peak. The amplitude for short 
distances from each peak decays very sharply, instead. 
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Fig. 3. The correlation function g(r) versus lnr in 2D for W = 0.1, 1, 3, 5, 7, 8, 9, 10, from bottom to 
top and system-size L = 121. All curves can be fitted by straight lines g(r) <~ 77 lnr + (3 and the power-law 
exponent r\ is plotted as a function of W in the inset. 

We have also obtained similar relations g(r) = 77 lnr + f3 for the cubic (3D) case by 
plotting g(r) vs. ln(r) for different values of disorder where g(r) still increases linearly 
with lnr. This means that the corresponding amplitudes show a power-law decay from the 
maximum, for any finite W, with the exponent 77 proportional to W. The unusual power-law 
localization of the mid-band state is rather surprising in 3D, since nearby states might be 
extended if the disorder is weak. In fact, 77 for both 2D, 3D, is found to increase almost 
linearly with W, shown in the insets of Fig. 3, 4, where we establish the very approximate 
relations 77 ~ W/3 in 2D and 77 ~ W/9 in 3D. 
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Fig. 4. The correlation function g(r) versus lnr for 3D with W = 0.1, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, from 
bottom to top and system-size L = 21. The curves can be fitted by straight lines g(r) <~ 77 lnr + [3 and the 
obtained r\ is plotted in the inset as a function of W. 

In summary, we have investigated the E = wave function for square and cubic lattices 
with off-diagonal disorder. This zero energy wave function in the adopted geometry exists 
for any disorder configuration and can be easily studied numerically. The wave function 
amplitude displays many sharp peaks randomly scattered in space and is characterized by 
the fractal dimension D 2 which shows a strong dependence on the strength of disorder. The 
amplitude from each peak falls off very rapidly for short ranges but the heights of the peaks 
decay slowly from the main maximum via the power-law r~ v , where 77 is almost linearly 
proportional to the strength of disorder W. 
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